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Abstract We consider the example from invariant theory concerning the conjuga- 
tion action of the general linear group on several copies of the n x n matrices, and 
examine a symmetric function which stably describes the Hilbert series for the in- 
variant ring with respect to the multigradation by degree. The terms of this Hilbert 
series may be described as a sum of squares of Littlewood-Richardson coefficients. 
A "principal specialization" of the gradation is then related to the Hilbert series of 
the K-invariant subspace in the GL„-harmonic polynomials (in the sense of Kostant), 
where K denotes a block diagonal embedding of a product of general linear groups. 
We also consider other specializations of this Hilbert series. 1 

To Nolan Wallach, who has influenced generations of mathematicians. 



1 Introduction 



We let F„ denote the finite dimensional irreducible representation of the general lin- 
ear group, GL„, (over C) whose highest weight is indexed by the integer partition 
A = (Ai > A2 > • • • > A„) (in standard coordinates). Given a finite sequence of in- 
teger partitions fi = (jlt^, • • • we will let be the multiplicity of in the 
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u (i) »H 

m-fold tensor product F„ ® • • • (g> F„ under the diagonal action of GL„. That is, 

denotes a (generalized) Littlewood-Richardson coefficient. 

We remark that the usual exposition of Littlewood-Richardson coefficients (see 
[Ful97,GW09,HL07,HL12,HTW05,Mac95]) concerns the case where m = 2. How- 
ever, by iterating the Littlewood-Richardson rule (or its equivalents) one obtains 
several effective combinatorial interpretations of our c^. 

The subject of this exposition concerns some interpretations of the positive in- 
teger £ (cjj) where the sum is over certain finite subsets of non-negative integer 
partitions. We believe that such sums have under-appreciated combinatorial signifi- 
cance. For example, one immediately observes the very simple specialization to the 
case where jit^) = (1) for all j = 1,- • • ,m, in which case the sum of squares is ml, 
which may be viewed as a consequence of Schur-Weyl duality. More generally, if 
V = (Vi > • • • > V,„ > 0) is a partition and n^> = (v 7 ), then is equal to the Kostka 
number Kx v (i.e. the multiplicity of the weight indexed by v in F^). 

Our motivation for considering these numbers comes from invariant theory. On 
one hand, we consider the conjugation action of the general linear group on sev- 
eral copies of the n x n matrices. On the other hand, we consider the K-conjugation 
action on one copy of the nx n matrices, where K denotes a block diagonal em- 
bedding of a product of general linear groups. These problems are related, and have 
been studied extensively. We make no attempt to survey the literature, but recom- 
mend [Dre07]. 

Central to this work is the notion of a Hilbert series. Let V be a graded vector 
space. That is, V = ©J = q Vd where is a finite dimensional subspace. The Hilbert 
series, £J = o(dimVrf)# d , formally records the dimensions of the graded components. 
Here q is an indeterminate. We also consider multivariate generalizations corre- 
sponding to situations where V is graded by a cone in a lattice. 

In our setting, V will be a space of invariant polynomial functions on m copies of 
the n x n matrices. For fixed values of the parameters, the Hilbert series is the Taylor 
expansion of a rational function around zero. When these parameters are small, one 
can expect to write down the numerator and denominator explicitly. These polyno- 
mials encode structural information about the invariants. However, as the size of 
the matrix becomes large, these rational functions are difficult to compute explic- 
itly. This motivates re-organizing the data by studying the coefficients of the Hilbert 
series of a fixed degree as the size of the matrix goes to infinity. The limit exists. 
The formal series recording this information will be referred to as the stable Hilbert 
series. 

Certain sums of squares of Littlewood-Richardson coefficients describe the co- 
efficients of the (stable) Hilbert series for the invariant algebra in each case. These 
Hilbert series, stably, may be expressed as a product. Furthermore, a "principal spe- 
cialization" of this product is then related to the Hilbert series of the K-invariant 
subspace in the GL„-harmonic 2 polynomials. 



2 Harmonic in the sense of Kostant (see [Kos63]), which generalizes the usual notion of a harmonic 
polynomial. 
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Unless otherwise stated, we will only need notation for representations with poly- 
nomial matrix coefficients, which are indexed by partitions with non-negative inte- 
ger components. The sum of the parts of a partition v will be called the size (de- 
noted |v|), while the number of parts will be called the length denoted £(v). As 
usual, we will also write A h d to mean |A| = d. Furthermore, we also adapt the 
(non-standard) notation that if fi = (jU^,- • ■ is a finite sequence of partitions, 

we set |/x| = YIj=\ IM^'I' an d write /x h d to mean = d. From a combinatorial 
point of view, the results involve a specialization of the following 

Theorem 1 (Main Formula). Let h,t2,h,--- denote a countably infinite set of in- 
determinates. We have 

°° 1 / \ 2 

fll-(rf + ^ + -) = ??W *" 

where the outer sum is over all partitions X and the inner sum is over all finite 
sequences of partitions // = (^ l \^ 2 \^ 3 \ ■ ■ •) with /.'' = '4^ '4^ ' ' ' '• 

Proof. See Section 7. □ 

As an application of the main formula, we turn to the space, Jt?(Ql n ), of GL„- 
harmonic polynomial functions on the adjoint representation (with its natural gra- 
dation) by polynomial degree, Jf(gl n ) = ®^ = q J^° d (gl n ) ■ The group GL„ acts on 
J^(gl„). Note that the constant functions are the only GL n -invariant harmonic func- 
tions. However, if K is a reductive subgroup of GL„, the space of K-invariant func- 
tions is much larger. Consider the example when the group K is the block diagonally 

embedded copy of GL ni x • • • x GL„ m in GL„, with n\ H Yn m =n. We will denote 

this group by K(n) where n — (m, • • ■ ,n m ). The purpose of this paper is to relate 
the dimension of the K(n) -invariant polynomials in J^f" i (gl n ) to a sum of squares 
of Littlewood-Richardson coefficients. See Theorem 2 for the precise statement. 

We consider this question because the related algebraic combinatorics are par- 
ticularly elegant, and hence have expository value in connecting harmonic analysis 
with algebraic combinatorics. However, this example is the tip of an iceberg. Indeed, 
one can replace GL„ with any algebraic group G (with g = Lie{G)) and K(n) with 
any subgroup of G. This area of investigation is wide open and well motivated as an 
examination of the special symmetries of harmonic polynomials. 

In Section 2, we describe a general "answer" to this question when K is a sym- 
metric subgroup of a reductive group G. This answer is not as explicit as would 
be desired, but applies to any symmetric pair (G,K). The remainder of the paper 
is related to the G = GL„ example with K = K(n). Note that when m = 2 (i.e. 
n = («i ,712)) the pair (G,K) is symmetric, but for m > 2 is not. We remark that in the 
m = 2 case, the results presented here were first described in [Wil07]. Our present 
discussion amounts to a generalization to m > 2. 

After setting up appropriate notation in Section 3 we provide an interpretation 
for a description of the Hilbert series of the K(n)-invariants in the GL„-harmonic 
polynomials on gl n in Sections 4 and 5. Chief among these involves sums of squares 
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of Littlewood-Richardson coefficients. We recall other combinatorial interpretations 
in Section 6. These interpretations involve counting the conjugacy classes in the 
general linear group over a finite field. 

Acknowledgements The first author wishes to thank Marquette University for sup- 
port during the preparation of this article. The second author wishes to thank the 
National Security Agency for support. We also would like to thank both Lindsey 
Mathewson and the referee for pointing out several references that improved the 
exposition. 



2 The case of a symmetric pair 

Let G denote a connected reductive linear algebraic group over the complex num- 
bers and let g be its complex Lie algebra. We have g = $(q) ®Q S s, where 3(0) denotes 
the center of g, while g ss = [g,g] denotes the semisimple part of g. A celebrated re- 
sult of Kostant (see [Kos63]) is that the polynomial functions on g, denoted C[g], are 
a free module over the invariant subalgebra, C[g] G , under the adjoint action. Choose 
a Cartan subalgebra \) of g, and let <P and W denote the corresponding set of roots 
and Weyl group, respectively. Choose a set of positive roots <P + , and let <P = <P + 
denote the negative roots. Set p = \ Y*ae<P + ct. For w eW, let l(w) denote the num- 
ber of positive roots sent to negative roots by w. Fix an indeterminate t . There exist 

positive integers e\ < ei < ■ ■ ■ < e r such that T.wew 1 '^ = 11;= 1 ^rH~ wnere r i s 
the rank of g ss . A consequence of the Chevalley restriction theorem ([Che55]) is 
that C[g] G is freely generated, as a commutative ring, by dim^g) polynomials of 
degree 1, and r polynomials of degree e\,-- ■ ,e r . These polynomials are the basic 
invariants, while e\, ■ ■ ■ , e r are called the exponents of G. 



2.1 Harmonic polynomials 

We define the harmonic polynomials on g by 

JT = {/ e C[g]\A(f) = for all A e D[ fl ] G } 

where D[g] G is the space of constant coefficient G-invariant differential operators 
on g. In [Kos63], it is shown that as a G-representation Jf° g is equivalent to the 
G-representation algebraically induced from the trivial representation of a maximal 
algebraic torus, T, in G. Thus, by Frobenius reciprocity the irreducible rational rep- 
resentations of G occur with multiplicity equal to the dimension of their zero weight 
space. Moreover, as a representation of G, the harmonic polynomials are equivalent 
to the regular functions on the nilpotent cone in g. 



Sums of squares of Littlewood-Richardson coefficients 



5 



The harmonic polynomials inherit a gradation by degree from C[g]. Set J^P* = 
n C[g]rf. Thus, = 0J = o We next consider the distribution of the mul- 
tiplicity of an irreducible G-representation among the graded components of J^. A 
solution to this problem was originally due to Hesselink [Hes80], which we recall 
next. 

Let fp t : h* — > N denote 3 Lusztig's q-analog of Kostant's partition function. That 
is fp t is defined by the equation, 

ae<f+ le 5ee(o,f)) 

where g(g, h) C h* denotes the lattice defined by the integer span of the roots. As 
usual, e^= denotes the corresponding character of T, with Lie(T) = h. As usual, we 
set^) = 0for^e( fl ,h). 

Let P(q) denote the integral weights corresponding to the pair (g, h). The dom- 
inant integral weights corresponding to <P + will be denoted P+(q). Let L{X) de- 
note the (finite dimensional) irreducible highest weight representation with highest 
weight X € P+(g)- The multiplicity of L(X) in the degree d harmonic polynomials, 
^ d {Q), will be denoted by nid(X). Set mx (t) = Y%=Q m d{X)t d . Hesselink's theorem 
asserts that 

m x {t)= I(-l)' w AWHp)-p). 

wew 

See [WWOO] for a generalization of this result. 

We remark that the above formula is very difficult to implement in practice. This 
is in part due to the fact that the order of W grows exponentially with the Lie algebra 
rank. Thankfully, only a small number of terms actually contribute to the overall 
multiplicity. See [Harll] for a very interesting special case where the number of 
contributed terms is shown to be a Fibonacci number. 



2.2 The K-spherical representations of G 

Let (G,K) be a symmetric pair. That is, G is a connected reductive linear algebraic 
group over C and K = {g e G\9(g) — g}, where 9 is a regular automorphism of G 
of order two. Since K will necessarily be reductive the quotient, G/K is an affine 
variety and the C-algebra of regular function C[G/K] is multiplicity free as a rep- 
resentation of G. This fact follows from the (complexified) Iwasawa decomposition 
of G. Put another way, there exists S C P+(q) such that for all X e P+(q) we have 



3 As always, N = {0, 1,2,3, • • •}» the non-negative integers. 
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Note that the subset S may be read off of the data encoded in the Satake dia- 
gram associated to the pair (G,K). The above fact implies that the Hilbert series 
H t (G,K) = Z7 =0 h d t d with h d = dim^ £/ (g) K has the following formal expression 

H t (G,K)= £ (-l)'M £ft(w(Hp)-p) . 

wew \ies J 

This formula seems rather encouraging. Unfortunately, the inner sum is very diffi- 
cult to put into a closed form for general w G W. This is, in part, a reflection of the 
fact that the values of p t cannot be determined from a "closed form" expression. 
However, note that w(X + p) — p often falls outside of the support of p t , and there- 
fore it may be possible to obtain explicit results along these lines. Moreover, the 
point of this exposition is to advertise that combinatorially elegant expressions may 
exist. At least this is the case for the pair (GL„ 1+ „ 2 ,GL ni x GL„ 2 ), as we shall see. 



3 Preliminaries 

We let q[ h denote the complex Lie algebra of n x n matrices with the usual bracket, 
[X,Y] = XY — YX, fovX,Y G qI„. Let Ejj denote the n x n matrix with a 1 in row 
;' and column j, and a everywhere else. The Cartan subalgebra will be chosen to 
be the span of {Eu\\ < i < n}. The dual basis in h* to (£n,£ , 22, ■ ■ ■ ,E nn ) will be 
denoted (£1, • • • ,£„). Choose the simple roots as usual, II = {e, — e (+ i 1 1 <;'<«}. 
Let (resp. <P + ) denote the roots (resp. positive roots). We will identify h with 
h* using the trace form, {H\,H2} = r Ir(H\H2} (for H\,H2 € h). The fundamental 
weights are a>; = Y!j=\ £j € h* for 1 < i < n- 1. We also set (On = £" =1 e 7 - e h*. Let 
P(GL„) = YJj=\ Z»y, and P+ (GL„) = Zco„ + ^l} Ncoj. 

From this point on, we will write (a\, ■■■,a n ) for Thus, we have X = 

(Ai, • • • ,X„) G P + (GL„) iff each A,- in an integer and X\ > ■ ■ ■ > X„. The finite di- 
mensional, irreducible representation of GL„ with highest weight X will be denoted 
(^> F «) wh ere 

K\ '■ GL„ ->• GL(F^). 

To simply notation will write Fjj - for (%,F^). 

Throughout this article, the representations of GL„ which we will consider have 
polynomial matrix coefficients. Thus the components of the highest weight X will 
be non-negative integers. Therefore, if X is a (non-negative integer) partition with at 
most I parts (£ < n), then the n-tuple, (Ai , • • • , Xg, 0, • • • , 0), corresponds to the highest 
weight of a finite dimensional irreducible representation of GL„ (with polynomial 
matrix coefficients). 

Under the diagonal action, GL„ acts on the c/-fold tensor product ® d( C n . Schur- 
Weyl duality (see [GW09] Chapter 9) asserts that the full commutant to the GL„- 
action is generated by the symmetric group action defined by permutation of factors. 
Consequently, one has a multiplicity free decomposition with respect to the joint 
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action of GL„ x Sd- Moreover, \fn>d then every irreducible representation of Sd 
occurs. The irreducible representation of Sd paired with will be denoted by v£. 
The full decomposition into irreducible GL n x ^-representation is 

X 

where the sum is over all non-negative integer partitions, X, of size d and length 
at most n. Note that when n — d, then the condition on £(X) is automatic. Thus, 
all irreducible representations of Vj occur. In this manner, the highest weights of 
GL„ -representations provide a parametrization of the ^-representations. 



3.1 Littlewood-Richardson coefficients 

Let d — [d\ , • • • , d m ) denote a tuple of positive integers with d = d\-\ h d m . Let 

Sd denote the subgroup of Sd consisting of permutations that stabilize the sets per- 
muting the first m\ indices, then the second ni2 indices etc. Clearly, we have 

S d =S di x---xS dm - 
The irreducible representations of Sd are of the form 

m(1) iiW 

where /x' 7 ' is a partition of size dj. It is well known that if an irreducible represen- 
tation, Vj, of Sd is restricted to Sd then the multiplicity of V(fi) in vf is given by 
the Littlewood-Richardson coefficient c^. This fact is a consequence of Schur-Weyl 
duality. 



4 Invariant polynomials on matrices 

A permutation of {1,2, • • • ,m} may be written as a product of disjoint cycles. This 
result is fundamental to combinatorial properties of the symmetric group, S m . Keep- 
ing this elementary fact in mind, let X — (Xi,X2,--- ,X m ) be a list of complex nxn 
matrices. Let Tr denote the trace of a matrix and define 

Tr a (X) = Tr(X (i)X (i) • • X m) Tr(X {2) X (2) • - X _ (2) )- ■ -Tv(X {k) X {k) ■ ■ X {k) ) 

(1) 

where a = (a^a^ ■ ■ ■ Om}){of^ ■ ■ ■ Om}) ■ ■ ■ (a^a^ ■ ■ ■ Om k ) is a permuta- 
tion of m = £m, written as a product of k disjoint cycles. Observe that the cycles 
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of (7 may be permuted and rotated without changing the permutation. In turn, the 
function Tr CT displays these same symmetries. 

Chief among the significance of Tr CT is the fact that they are invariant under the 
conjugation action g ■ (X u - ■ ■ ,X m ) = {gXig~ 1 , • • • ,gX m g~ 1 ) , where g € GL„. By set- 
ting some of the components of (X\ ,X m ) equal one defines a polynomial of equal 
degree but on fewer than m copies of M„. Intuitive, this fact may be described by 
allowing equalities in the components of the cycles of a. That is to say (formally) 
we consider a up to conjugation by Levi subgroup of S m . 

In [Pro76,Pro76b], C. Procesi described these generators for the algebra of GL„- 
invariant polynomials, denoted 4 C[V] GL ", on V = M™, and provided a proof that 
these polynomials span the invariants. Hilbert tells us that the ring of invariants must 
be finitely generated. Thus, there must necessarily be algebraic relations among this 
(infinite) set of generators. In light of Procesi's work, these generators and rela- 
tions are understood. However, recovering the Hilbert series from these data is not 
automatic. 

In order to precisely quantify the failure of the Tr CT being independent, we intro- 
duce the formal power series A n (t ) = A„ (t i , t 2 , • • • , t m ) = £ a n (d) t d with 5 the coeffi- 
cient defined as a„(d) = dimC[M™]j L " where C[M™] d = C[M n © • • • ®M n \ dl ... 4m) 
denotes the homogeneous polynomials of degree dt on the ;'th copy of M n . The mul- 
tivariate series, A n {t), is called the Hilbert Series of the invariant ring. Except in 
some simple cases, a closed form for A n (t) is not known. Part of the present expo- 
sition is to point out the rather simple fact that a n (d) may be expressed in terms of 
the squares of Littlewood-Richardson coefficients. Furthermore, we prove 

Proposition 1. For any natural numbers m and d= (d\ , • • • ,d m ) 

the limit lim„^oo a n (d) exists. If we call the limiting value a{d) and set A(t) = 

Y<a{d)t d then 

~ 1 

A(t) = Hi-(^+--+4)' 

Proof. The polynomial functions on M n are multiplicity free under the action of 
GL„ x GL„ given by {gi,g 2 )f(X) = f(gi l Xg 2 ) for g ug2 G GL„, X e M n and 
/ € C[M„}. The decomposition is a "Peter- Weyl" type, C[M„] = ©(F^)* ® F^, 
where the sum is over all non-negative integer partitions X with l(X) <n. We have 
C[®? =1 M n ] S* ®™ =1 C[M„]. Thus 




4 Here we let M„ denote the complex vector spaces ofnx/i matrices (with entries from C). 

5 We will use the notation d = (d\ , ■ ■ ■ , d m ) , t d = ff 1 • • • f*» . 
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with respect to the GL„ x GL„ action on the diagonal of M'". Note that in multide- 
gree (d\, ■ ■ ■ ,d m ) the sum is over all fi with \fi^ \ — dj. If we restrict to the subgroup 
{(SiS) '■ 8 <= GL„} of GL„ x GL„, we obtain an invariant exactly when a = j3. The 
dimension of the GL„ -invariants in the degree d homogeneous polynomials on M™ 

is therefore £ (c^) where the sum is over all 11 = (ju (1) ,- • • ,ju (m) ) and X h m with 
length at most n. The degree d component decomposes into multi-degree compo- 
nents (di , ■ ■ ■ ,d m ) with d — Y^dj. 

If d >n, then the condition that £(X) < n is automatic, and this fact implies that 
if 7^ 0, then £{p.^) < n for all j. Thus, the dimension of the degree d invariants 

in C[M™] is £ (c^) 2 where the sum is over all partitions of size d. If we specialize 
the main formula so that t j = for j > m then the sums of squares of Littlewood- 
Richardson coefficients agree with A(t). □ 

For our present purposes, we will specialize the multigradation on the invariants 
in C[M™] to one that is more coarse. From this process, we can relate the stabilized 
HUbert series of the invariants in Jf?[M n ]. This specialization will be the subject of 
the next section. 



We now turn to another problem that, we shall see, is surprisingly related. Con- 
sider the nxn complex matrix 



X 



X(l,l) X(l,2) ■■■ X(l,m) 
X(2,l) X(2,2) ••• X(2,m) 



_X(m,l) X(m,2) ■■■X(m,m) 
where X(i,j) is an m x n ; complex matrix with n = Y J rij. Define 

Tr CT (X) = f\Tr(x(a[-^4 j) )X(4 J \^ j) )X(^ J \ai j) )---X(a^ . 
j=i v ' 

Let K(n) denote the block diagonal subgroup of GL„ of the form 

GL ni 



K(n) 



GL «„, 



The group K(n) acts on M n by restricting the adjoint action of GL„. The K(n)- 
invariant subring of C[M„] (denoted C[M„} K{n ^ is spanned by Tr a (X). For small 
values of the parameter space, these expressions are far from linearly independent 
(as in the last example). Formally, one cannot help but notice the symbolic map 
Tr a i ^ Tr<j. We will try next to make a precise statement along these lines. 
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Define C[M n ]</ to be the homogeneous degree d polynomials on M n , and let 
C[M n ]f n) denote the K(n)-invariant subspace. Set a^ n \d) = dimC[M n ]f n \ and 
A( n \t) = T.°d=o<^ n \d) t d . Analogous to Proposition 1, we have 

Proposition 2. For any non-negative integer d, the limit 
lim lim ••• lim a^ nu '"' nm \d) 

n\— 5-oo«2— 5-°° n m — 

exists. Denote the limiting value a(d) and set A{t) = Y.d=0 a (d)t d - We have 

Proof. We begin with the GL„ x GL„ -decomposition C[M„] = (F*)* ®F^ with 
respect to the action in the proof of Proposition 1 . The irreducible GL„ -representation 
is reducible upon restriction to K(n). The decomposition is given in terms of 
Littlewood-Richardson coefficients 

Therefore, as a K(n) x K(n) -representation we have 

C[M„] rf = £ (p^v) (®7=iFjf)*® (®7=i F »f ) 

where the sum is over all A with |A| = J, £(A) < « and £(n^),£(v^) < tij. If we 
restrict to the diagonally embedded K(n)-subgroup, we obtain an invariant exactly 

when (a = v. Thus, a^ n \d) = £ (c^) , with the appropriate restrictions on X and 

A 4 - 

If all rij > d then the condition on the lengths of partitions disappears, and we 
may sum over all A h d. The result follows by specializing the main formula by 
setting tj — t for 1 < j <m and tj = for j > m. □ 

Although we will not need it for our present purposes, it is worth pointing out 
that the algebra C[M n ] has a natural N m gradation defined by the action of the center 
of K(n). This multigradation refines the gradation by degree. The limiting multi- 
graded Hilbert series is the same as that of Proposition 1 . Upon specialization to 
t\ = t2 = ■ ■ ■ = t m = t we obtain the usual gradation by degree (in both situations). 
The advantage of considering the more refined gradation is that one can consider 
the direct limit as m and n go to infinity. This will be relevant in the next section. 
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5 Harmonic polynomials on matrices 

A specific goal of this article is to understand the dimension of the space of degree d 
homogeneous K(n)-invariant harmonic polynomials on M n . With this fact in mind, 
we observe the following specialization of the product in the main formula. Let 
tj = t>. Then we obtain: 

oo j ooj oo l_fk 

n t^tt^tt^th = n tt^z: = n tt^i • 

For a sequence fi = {fi^ l \^ 2 \ ■ ■ ■) set gr(/j,) — £°° =1 j\n^\. The equation in the 
main formula becomes 

n^=Li(4) 2 ^ ) - 

k=\ 1 zr X fj. 

The notation gr is used for the word grade. We explain this choice next. Let C [M„ ; d] 
denote the polynomials functions on the nxn complex matrices together with the 
gradation defined by d times the usual degree. That is, <C[M n ;di]^ denotes the 
degree di homogeneous polynomials, but regarded as the d\di graded component 
in C\M n \d\]. We consider the N-graded complex algebra srf n defined 6 as 

s/n - C[M„;1]®C[M„;2]®C[M„;3]®--- 

oo 

5=0 

where srf n \8\ is the graded 8 € N component (with the usual grade defined on 
a tensor product of algebras). The group GL„ acts on each C[M n ;d] by the ad- 
joint action, and respects the grade. Under the diagonal action on the tensors, the 
GL„-invariants, ^ GL ", have A„(t) as the Hilbert series when t is specialized to 
(fi,f2,?3, • • •) = (f ,f 2 ,f 3 , • • •). This GL„-action respects the 8 component in the gra- 
dation. That is to say that the Hilbert series is 

A„(f,f 2 ,f 3 ,-..)= £dimK[5]) GL "f 5 . 

5=0 

For d — 0, • • • ,n the coefficient of t d in A n (t,t 2 , • • •) is the same as the coefficient of 
t d in A(t ,f 2 , • • •). Summarizing, we can say that stably, as n — >• °°, the Hilbert series 

of< GL "isnr=ii^- 

We turn now to the GL„-harmonic polynomials on M„ together with its usual 
gradation by degree. Kostant's theorem [Kos63] tells us that 



6 One might consider this as a "Bosonic Fock space." Precisely, we mean the infinite tensor product 
in which each element is a finite sum of tensors involving only finitely many of the factors. 
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C[M n ] ~ C[M„] GL " <g> (Af„). 



As before let K = K(«i,«2) denote the copy of GL ni x GL„ 2 (symmetrically) em- 
bedded in GL„ |+ „ 2 . Passing to the K(m,tt2) -invariant subspaces we obtain: 

C[M„f = C[M n f L "<E)Jf(M n ) K . 

The Hilbert series of C[M„] GL " is well known to be IlLi rrr> while the Hilbert 
series C[M„] K is A„(t,t). These facts imply that the dimension of the degree d 
homogenous K-invariant harmonic polynomials is the coefficient of t d in F n {t) = 
A„(/,()[]" =1 (l - t J ). We have that the coefficient of f d in F„(t) for 
<i = 0, • • ■ ,min(m,n2) agrees with the coefficient of t d in 

^ 00 00 1 A 

k=\ k=l 1 

Again summarizing, we say that stably, as «i , «2 — >■ 00 , 

the Hilbert series of Jf? (M„ 1+ „ 2 ) K ( Ml '" 2 ) is the same as ,< GL " as n,m,n 2 -> °°. That 
is to say, for fixed 8 we have 

lim dimJf(M n . +n? )f" [ '" 2) = lim dim^/„[5] GL ". 



Observe that this procedure generalizes. Let m > 2. Analogous to before, let 
C[M™;<f] denote the polynomial function on M™ = M n © • • • ®M n (m-copies) to- 
gether with the N-gradation defined such that C[M™;d\]ci 2 consists of the degree 
d2 homogeneous polynomials but regarded as being the (ficfe-th graded component. 
Note that the 5-th component in the grade is (0) if 5 is not a multiple of d\. 

Let srf™ be the N-graded algebra defined as 

<" = C[M™- l -\]®C[M™- l -2] ®C[AC _1 ;3] <8> ■ ■ ■ 

Since stf™ is a tensor product of N-graded algebras, it has the structure of an N- 
graded algebra. As before let srf™\?>\ is the 5-th graded component. The group GL„ 
acts on each C[M™;d] by the adjoint action, and respects the grade. 

Next, set t\ = ti = ■ ■ ■ = f m _ 1 = t , then t m = t m+ \ =■■■ = f 2m -2 = t 2 , then ti m - 1 = 
• • • = ?3m-3 = t 3 , and so on. The result of this procedure is 

Theorem 2. For all n = (m, • • •,«,„) e Z™ and rfeN, Zer 

/ , \ K(n) 

^(n) = dim(^(M„)J . 

I/zen for fixed d, the limit 



lim ••• lim h d (n u - ■ ■ ,n m ) 
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exists. Let the limiting value be h^. Then 

h d = limdim« m [d]) GL \ 

n— 5.00 

Proof. After the specialization we obtain: 

n i_(f*+... + f* + f2i + ... +f 2t + ... ) - n 1 _( OT _ 1 )( f t +r 2* + ...) 



m-i copies m-i copies 



00 J 



1 



. . t k -(m-\)t k 



f,_i 1 



l-t k 



k=l 



mt 



k ■ 



The significance of this calculation is that it allows for another interpretation of 
sums of Littlewood-Richardson coefficients. The rest of the proof is identical to the 
m = 2 case in the preceding discussion. □ 



5.1 A bigraded algebra and a specialization of the main formula 

As before, the group GL„ acts on M n by conjugation, and then in turn acts diagonally 
on M n ®M n . That is, given (X,Y) € M„ ®M n , and g E GL„, we have g ■ (X,Y) = 
(gXg~ 1 , gYg~ 1 ) . We then obtain an action on C [M„ © M n ] . 

We have already observed that the algebra C[M„ S)M n ] is bigraded. That is, let 
C[M n ®M n ](i,j) denote the homogenous polynomial functions on M n (BM n of de- 
gree i in the first copy of M n and degree j in the second copy of M n . Let a and b be 
positive integers. We set 

C[M n M„; (a,b)] (ai, bj) = C[M„ @M n ] (i, j) 

with the other components zero. Put another way, C[M„ ffiM„; (a,b)] is the algebra 
of polynomial function on M n ®M„ together with the bi-gradation defined by (a,b) 
times the usual degree. As before we consider the infinite tensor product 

00 00 

&n = C \M n M n ; ( a , b ) ] . 

a=16=l 

Next, note the following, obvious, identity: 
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n- 

k=\ l 



(l-^)(l-y*) 



n 

k=\ 



{l-x k )(l~y k ) 
\-{x k +y k ) 



We observe that this is a specialization of the product in the main formula. Specifi- 
cally, let q'tf = z s where z s is given as the entry in the following table. 





q q 2 q 3 q 5 q 6 ■■■ 


t 


Z\ Z2 Z4 Zl Zll ••• 


t 2 


Z3 zs z% zn zn ■■■ 


r 3 


Z6 Z<) Z13 Zl8 Z2A 


f 4 


ZlO ^14 Z19 Z25 Z32 


f 5 


Zl5 ^20 Z26 Z41 ' ' ' 



Then, 



£=1 1 



(4+4- 



-A- 



r nr ^ 



n; 



k=l 



J 



= n oo (l-/)(l-r^) 
i «=l l-(^+f*) ' 

In this way, the GL„ -invariants in the harmonic polynomials on M n Q)M„ may be 
related to the multigraded algebra structure as in the singly graded case. 

This identity becomes significantly more complicated when generalized to har- 
monic polynomials on more than two copies of the matrices. This fact will be the 
subject of future work. 



6 Some combinatorics related to finite fields 

In this section, we collect remarks of a combinatorial nature that provide a more 
concrete understanding of the sum of squares that we consider in this paper. 
From elementary combinatorics one knows that the infinite product 

k=l 1 < l t n=Q£=0 



where p n j is the number of partitions of n with exactly I parts. When q is specialized 
to a positive integer the coefficients of this series in t has many interpretations. 
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6.1 The symmetric group 

Fix a positive integer d, and a tuple of positive integers d = {d\,d%, - ■ ■ ,d m ) with 

d\-\ h d m = d. As before, let 5d denote the subgroup of S d isomorphic to 5^ x 

• • • x Sd m embedded by letting the i ,h factor permute the set /, where {J\, ■ ■ ■ ,J m } 
is the partition of {1,- • • ,d} into the m contiguous intervals with = d[. That is, 
J\ = {l,2,---,di}, h = + -,di+d 2 }, etc. 

The group 5d acts on 5^ by conjugation. Let the set of orbits be denoted by ff(d). 
We have 

Proposition 3. For any d = {d\ , ■ ■ ■ ,d m ), 

where the inner sum is over all tuples of partitions with jX^ h dj. 
Proof. We begin with the Peter- Weyl type decomposition of C [Sd] 

C[5,/]=© (v d x )*®V d \ 

Ah d 

We then recall that Littlewood-Richardson coefficients describe the branching rule 
from Sd to Sd- 

Combining the above decompositions, we obtain the result from Schur's Lemma. 
□ 

It is an elementary fact that the Sj-conjugacy class of a permutation a in Sd 
is determined by the lengths of the disjoint cycles of a. A slightly more general 
statement is that the Sd-conjugacy class of a € S d is determined by a union of cycles 
in which each element of a cycle is "colored" by colors corresponding to J\, ■ ■ ■ ,J m . 
It is not difficult to write down a proof of this fact, but we omit it here for space 
considerations. 

If one fixes d, then the sum over all <f -compositions, £d|^(d)|, may be inter- 
preted as the number of J-bead unions of necklaces with each bead colored by m 
colors. For example ifm — 2, and d = 4, the resulting set may be depicted as: 
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9999 
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9999 
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<0 o 
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9999 
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oo 
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9 9 
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oo 


f>oO 


9999 



A single fe-bead necklace has Z/^Z-symmetry. If the beads of this necklace are 
colored with m colors then the resulting colored necklace may have smaller group 
of symmetries. Using Polya enumeration (i.e. "Burnside's Lemma"), one can count 
such necklaces by the formula 

Nk(m) = \ Y<${ r ) m ~ r - 

K r\k 

where <j) denotes the Euler totient function 7 . The theory forming the underpinnings 
of the above formula may be put in a larger context of cycle index polynomials. We 
refer the reader to Doron Zeilberger's survey article (IV. 18 of [GBGL08]). 

The generating function for disjoint unions of such necklaces can be given by the 
product 



7 That is, 0(r) is the number of positive integers relatively prime to r. 
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» / 1 V 



That is to say that the number of c/-bead necklaces, counted up to cyclic symmetry, 
is equal to the coefficient of t d in 77(f). The main formula specializes to 



Mt)= E (<i) 2 tw, (2) 

where the sum is over partitions X and m-tuples of partitions, [i = (jU^, • • • 

Equation 2 begs for an (explicit) bijective proof which is no doubt obtained by 
merging both the Littlewood-Richardson rule and the Robinson-Schensted-Knuth 
bijection. It is likely that more than one "natural" bijection exists. 



6.2 The general linear group over a finite field 

Let p denote a prime number, and v € Z + . Set q = p v . Let GL m (q) denote the 
general linear group of invertible mx m matrices over the field with q elements. 
The set, %„(q), of conjugacy classes of GL,„(q) has a cardinality of note in that the 
infinite products has the following expansion: 

00 1 — t k 

k=l 1 t l t m=0 

see [BFH86, Sto88], and the references within. 

Note that from Equation 2 we obtain a new formula for the number of conjugacy 
classes of GL m (q). Namely, 



na-<'>Wn(r^) 



N k (q)-l 



Stong reference 



7 Proof of the main formula 



Before proceeding, we require some more notation. Given a partition X = (Ai > 
X2 > ■■■) with vi ones, and V2 twos, etc., we will call the sequence (vi, V2, • • •) the 
type vector of X. Note that £(X) — £v;, while the size of X ism = \X \ = = £zV;. 
As is standard, we set 
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It is elementary that the cardinality of a conjugacy class of a permutation with cycle 
type X is — . (Equivalently, the centralizer subgroup has order z\.) 

We now prove the main formula. The product on the left side (LS) may be ex- 
panded using the sum of a geometric series, the multinomial theorem, and then the 
sum-product formula, as follows 



co i oo oo 

" = n 1 _ (tt+<t+ , t+ ,„ ) = n£(^+-r 

k=l 1 (,r 1 +r 2 +r 3+ ) k=lu=0 



k=lu=0 «i+«2+-=« 1 2- 

ny («1 +M2-I )- t kui ku 2 
, U\\U 2 \--- 1 2 

k=i «i,«2,— 1 z 

00 fi/') V ■ (0 • (0 

M „(0,„(0, 1 2 "' 

(i) i=l Mi !M, ! ■ ■ ■ 

« ; . ,- 

(with all sequences having finite support.) We will introduce another sequence, a = 
(fli,fl2, • • •) and extract the coefficient of t a — Y\tf' in the above formal expression 
to obtain 

ETT ( M l + K 2 ^ 
11 (i). (i) 



„(i).wv V ;„(')_„ ' = 1 ^1 ;t *2 : " 



The above is not such a complicated expression, although these formal manipu- 
lations may, at first, seem daunting. Observe that the sequence uj\u^\--- with 
£i iuf — «/ encodes a partition of aj with the number i occurring exactly times, 
while X)i>i is the length of the partition. We will call this partition jU^. That is 
to say, jU ^ has type vector {uj ^ , w?' ,•••)• The coefficient of i a may be rewritten as 
a sum over all double sequences such that for all j, j'm^ = a 7 of 

A (4'V4 + --0i _n (» ( i +4 + --0i' E " y) m 

Given a (finitely supported) sequence of partitions fi = fj,( 2 \ ■ ■ •), we denote 
the partition obtained from the (sorted) concatenation of all jx^ by \JjX^\ Thus, if 

jU^) has type vector (u^ , , • • •) then the number of z"s in Ufi^ is + u\ ' H ■ 

It therefore follows that the numerator of the righthand side of Equation 3 is zx when 
A = UjU^). The denominator can easily be seen to be z^yj. From this observation 
we obtain 

TS _y Z ^ U) ><"|>< 2) I>< 3 >I 

u_ Lrr 7,, Zl 2 3 
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7.1 An application of the Hall scalar product 

Let A n denote the S„-invariant polynomials (over C as usual) in the indeterminates 
x\ , ■ ■ ■ ,x n , let A [x] = lim^_ A„ denote the inverse limit. Thus, A [x] is the algebra of 
symmetric functions. For a non-negative integer partition X, we let sx (x) denote the 
Schur function. That is, for each n, sx (x) projects to the polynomial sx {x\ , ■ ■ ■ ,x n ) 
which, as a function on the diagonal matrices, coincides with the character of the 
GL„-irrep . The set {sx (x) : X h d} is a C-vector space basis of the homogeneous 
degree d symmetric functions. We will define a non-degenerate symmetric bilinear 
form (•, •) by declaring (s a (x),sp (x)) = 8 a p for all non-negative integer partitions 
a, j3. This form is the Hall scalar product. 

Given an integer m, define p m (x) = x™ +x™ H to denote the power sum sym- 
metric function. Given v h N, set p v {x) — Y\Pvj (x). We remark that the left side of 
the main formula is easily seen to be Y,vPv(t)- 

It is a consequence of Schur- Weyl duality that the coefficients of the Schur func- 
tion expansion 

x 

are the characters of the 5^-irrep indexed by X evaluated at any permutation with 
cycle type v. It is a standard exercise to see, from the orthogonality of the character 
table, that Pv(x) are an orthogonal basis of A [x], and (p v (x),p v (x)) = z v . 

We next consider another set of indeterminates, y = yi,y2, "- Set A[x,y] = 
A [x] <gi A [y]. The Hall scalar product extends to A [a?, y] in the standard way as 

(f(x) ® g(y),f'(x)g' (y)) = (f(x)j'(x)) (g(y),g'(y)) . 

where f(x),f'(x) <£ A[x] and g(y),g'(y) <E A[y\. (We then extend by linearity to 
all of A[x,y].) 

The character theoretic consequence of (GL„,GL,t)-Howe duality is the Cauchy 
identity: 

i 

i,j=\ 1 x 'y.i x 

(in the infinite sets of variables). In fact, for any pair of dual bases, ax,bx, with 
respect to the Hall scalar product, one has nij=i l-xy ~ ^X a x( x )bx(y). From 
this fact one obtains 

1 

n ~rz — = l J Px{x)Px{y)lzx- (5) 

ij=i 1 Myj x 
From our point of view, we will expand the following scalar product 
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in two different ways, corresponding to Equations 4 and 5. 

First by homogeneity of the Schur function and Cauchy's identity 

nE'c^rf = e n^«>wnv)Wi (1)| ^ (2,L - 

Since the multiplication of characters is the character of the tensor product of the 
corresponding representations, we have s a sp — L c ^ s / ™ tne x ( res P- V) variables. 
Expanding the above product gives 

(nnd^.nd^)=E(4) 2 <r i «r i ---. 

\ifc=i u=i 1 x, yj' k u=i 1 / m 

Secondly, the scalar product in (6) may be expressed as 

(ft L/y^W/y^yVv*) t |,4W 1 , (w)Aa )■ 

\*=1 n(k) x I 

We observe that nr=i P M W = PA where ^ = U^ =1 £i W - 



7.2 A remark from "Macdonald's book" 

The results presented in this paper emphasize describing the cardinality of an or- 
bit space by a sum of Littlewood-Richardson coefficients. It is important to note, 
however, that the main formula can be written simply as 

Em*)=e(e(^)VVV 3) --- 

v \ X J 

So from the point of view of [Mac95], one realizes that the main formula is sim- 
ply a way of expanding Y.Pv{x). With remark in mind, we recall the "standard" 
viewpoint. 

For a non-negative integer partition 5 = (5i > 52 > • • •) let x s = xf'x^x^ 
The monomial symmetric function, mg(x), is the sum over the orbit obtained by all 
permutations of the variables. The monomial symmetric functions are a basis for the 
algebra A . 

The question becomes obtaining an expansion of into monomial sym- 

metric functions. This question is answered immediately by observing the expansion 
of py into monomial symmetric functions, which can be found in [Mac95] on page 
102 of Chapter I, Section 6. For partitions y and 8 define L(y,8) by the expansion 



Pr( x ) =Y< L ys m 8( x )- 

5 
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We next provide a combinatorial description of Lyg- Let 7 denote a partition of 
length I. Given an integer valued function, /, defined on {1, 2, 3, • • • ,£}, set 

f(Y)i= E vj 

j:fU)=i 

for each ; > 1 . 

The sequence (/(t)i,/(7)2,/(7)3, ■ ■ ■) does not have to be weakly decreasing. 
For example, if y = (1,1,1) and /(l) = 1, f(2) = 4 and /(3) = 4 then /(y)j = 1, 
7(7)4 = 2 and /(y)t = for all k 7^ 1,4. However, often this sequences does define 
a partition. We have 

Proposition 4. Ly § is equal to the number of functions f such that f(y) = 8. 

Proof See [Mac95] Proposition I (6.9) □ 

From Proposition 4 and the main formula, we obtain 

Corollary 1. Given a partition 8, the cardinality of 

{f\for some parition 7, /(y) = 8 } 

is equal to 

EE(^) 2 

Am 

where the sum is over all (J, such that jjU^ | = Sjfor all j and \X \ — 
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